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ABSTRACT 

This  research  was  concerned  widi  a  practical  approach  to  the  design  of  an  improved  age- 
hardenable  aluminum  alloy  for  moderate  temperature  application.  The  process  involved 
extensive  empirical  research,  theoretical  snnulation  and  modeling,  calculated  phase  diagrams  and 
first  principle  calculations  and  analytical  techniques.  Microstructure  characterization  and 
analytical  transmission  electron  microscopy  were  employed  to  refine  the  development  of  the 
quaternary  Al-Cu-Mg-Ag  phase  diagram.  Differential  scanning  calorimetry  and  electron 
diffraction  were  implemented  to  investigate  the  microstructural  evolution  of  the  ternary  alloys 
and  to  validate  calculated  equilibrium  phase  boundaries.  Energy  dispersive  spectroscopy  was 
also  used  to  illustrate  the  effect  of  trace  additions  on  the  phase  boundaries  that  are  important  for 
the  thermal  stability  of  this  class  of  alloys.  Computer  simulation  and  modeling  were  used  to 
identify  the  microstructure  that  will  result  in  the  optimal  combination  of  mechanical  properties. 
Unavailable  parameters  were  calculated  from  first  principle  atomistic  modeling.  The  major 
objective  of  this  research  program  has  been  to  illustrate  the  use  of  modeling  and  simulation  in 
alloy  design  and  therefore  aid  in  tiie  early  insertion  of  new  high  performance  materials. 

INTRODUCTION  AND  TECHNICAL  BACKGROUND 

Due  to  economic  considerations,  alloy  design  and  development  can  no  longer  be 
performed  using  the  purely  empirical,  trial  and  error  approach  that  has  been  dominate  over  the 
past  100  years.  Modeling  and  simulation,  in  conjunction  with  experiments,  can  be  employed  to 
improve  the  efficiency  of  alloy  design,  optimizing  processing  and  manufacturing  operations. 
This  approach  can  greatly  reduce  the  time  associated  with  alloy  development  and  the  generation 
of  the  knowledge  base  required  for  insertion  of  new  materials  into  aerospace  systems.  This 
research  and  development  program  was  aimed  at  obtaining  a  high-strength,  thermally  stable 
aluminum-based  alloy  that  would  be  suitable  for  moderate  temperature  applications  up  to 
approximately  150°C,  that  could  be  used  for  the  fuselage  of  a  Mach  2.2-2.4  aircraft.  Ultimately, 
the  objective  of  fiiis  research  is  to  provide  an  example  of  streamlined  alloy  design  based  upon  the 
marriage  of  theoretical  metallurgical  principles,  calculated  phase  diagrams  and  experimental 
alloy  design  and  processing.  An  outline  of  the  approach  undertaken  for  this  program  is 
schematically  represented  in  Figure  1. 


Figure  1  -  Alloy  development  methodology 

The  first  step  in  alloy  design  is  to  select  a  system  that  offers  promise  for  obtaining  the 
desired  microstructure  and  properties.  After  an  extensive  literature  search  and  examination  of 
empirical  data,  we  selected  the  Al-Cu-Mg-(Ag)  system  for  our  base  system.  In  order  to  identify 
the  optimized  precipitate  structure  and  morphology,  our  research  team  developed  a  dislocation- 
slip  simulation  method  [1],  In  combination  with  a  quantitative  microstructure  characterization 
procedure,  we  have  been  able  to  evaluate  superposition  of  the  strengthening  effects  of 
precipitates  with  different  volume  fractions,  sizes,  shapes  (aspect  ratio)  and  orientations  with 
respect  to  the  Al-matrix.  This  has  allowed  us  to  find  the  optimized  precipitate-structure  for 
maximum  strength.  For  instance  in  an  Al-Li-Cu  alloy  where  the  main  strengthening  precipitates 
are  Ti-{I1I}  and  e'-{100}  plates,  a  structure  comparison  factor  can  be  introduced  to  describe 
the  superposition  effects. 
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where  the  coefficient  ^  accounts  for  the  effect  of  plate  orientation,  and  are  the  mean 
diameter  and  the  aspect  ratio  of  the  A-plates,  respectively,  and  Dg  and  y^  for  those  of  the  B- 
plates.  We  have  shown  that  co-existence  of  both  T,  {111}  and  0'  {100}  plates  in  certain 
combinations  has  the  highest  overall  strengthening  effect  for  a  given  Cu-composition  [1].  The 
knowledge  gained  through  the  simulation  has  been  verified  in  some  Al-alloys  [1,2].  Regarding 
the  Al-Cu-Mg-Ag-X  system,  have  been  concerned  with  Q  (AljCu),  a  plate-like  phase  with  a 
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{111}„  habit  plane,  GXAljCu)  plate-like  phase  with  a  {100}„  habit-plane,  and  S  (AljCuMg)  a 
rod-like  phase  oriented  along  the  <210>„  direction.  These  three  phases  may  co-exist  in  many  Al- 
Cu-Mg-Ag  alloys  under  conventional  processing  conditions.  The  rod-like  S  or  S'  phase  is  less 
effective  in  the  strengthening,  possibly  due  to  its  particular  orientation,  which  results  in  relatively 
small  obstacles  to  dislocation  motion  on  {111}„  slip  planes.  Our  dislocation-slip  simulations 
have  shown  that  for  optimal  strength,  the  coexistence  of  a  fine  uniform  distribution  of  6  and  Q 
plates  in  certain  combinations  is  desired. 

The  Al-Cu-Mg-(Ag)  System:  We  selected  this  system  after  an  examination  of  the  literature  and 
available  empirical  data.  The  strengthening  precipitates,  especially  the  {100}„  AhCu  phase 
(0V0),  in  the  Al-Cu-Mg-X  2xxx  alloy  system  are  more  thermally  stable  at  moderate  temperatures 
than  those  present  in  the  Al-Zn-Mg-X  7xxx  or  the  Al-Si-Mg  6xxx  alloy  systems.  The  addition 
of  small  amounts  of  Ag,  on  the  order  of  0.4  '^/o,  to  Al-Cu-Mg-X  alloys  with  high  Cu:Mg  ratios 
(e.g.  10:1  '^/o)  has  been  shown  to  significantly  enhance  the  age  hardening  response  over  that  of  a 
similar  ternary  alloy  (c.f.  [3]).  The  marked  gain  in  strength  is  due  to  the  formation  of  a  uniform 
dispersion  of  a  hexagonal  plate  precipitate  with  a{lll}„  habit  plane,  designated  Q.  Several 
detailed  studies  have  indicated  that  the  internal  crystal  structure  of  die  £2  plates  is  quite  similar  to 
that  of  equilibrium  0  ([4-8]).  As  will  be  descriW  later,  we  have  assessed  the  validity  of  this 
assumption. 

It  has  recently  been  shown  that  Ag  facilitates  Ag-Mg  clusters  that  in  turn  act  as 
nucleation  sites  for  the  Q  phase  [9,10].  Detailed  atom  probe  and  high  resolution  Z-contrast 
microscopy  have  revealed  that  Ag  and  Mg  are  present  at  the  coherent  interfaces  between  the  Q 
plate  and  the  enveloping  aluminum  matrix  [9-11].  Two  prominent  factors  that  jeopardize  the  fine 
uniform  distribution  of  Q  plates  in  the  quaternary  alloys  have  also  been  identified.  The  long¬ 
term  stability  of  the  Q  plate  is  sensitive  to  the  presence  of  S-phase  precipitation  (Al2CuMg).  If 
the  alloy  composition  lies  within  the  equilibrium  (a  +  0  +  S)  phase  field  the  existence  of  S-phase 
will  create  chemical  potential  (p)  gradients  in  Mg.  Essentially  the  presence  of  S-phase  limits  the 
Mg  available  for  the  coherent  interfaces  of  the  Q  plates.  Similarly,  several  studies  have 
illustrated  Ag  segregation  to  S-phase  in  this  class  of  quaternary  alloys  [10,12].  The  presence  of 
S-phase  is  deleterious  to  the  stability  and  coarsening  resistance  of  Q  plates  due  to  the  fact  that 
the  Q  plates  may  only  exist  with  Mg  and/or  Ag  at  the  coherent  interfaces.  The  first  goal  was  to 
design  a  quaternary  alloy  that  lies  within  the  (a  +  0)  phase  field  for  a  wide  range  of  desirable 
temperatures.  Therefore,  we  must  have  a  thorough  understanding  of  phase  space  for  both  the  Al- 
Cu-Mg  ternary  and  the  Al-Cu-Mg-Ag  quaternary  Al-rich  alloys.  Once  within  the  (a  +  0)  phase 
field  the  driving  force  for  coarsening  of  the  Q  plates  must  be  quelled  in  order  to  maintain  a  fine 
uniform  and  thermally  stable  distribution. 

Phase  Diagram  and  Thermodynamic  Modeling:  The  CALPHAD  method  (CALculation  of 
PHAse  Diagrams)  is  an  approach  to  establishing  equilibria  among  phases,  through 
thermodynamic  modeling  of  the  individual  phases  based  on  the  idea  of  phase  competition  in  a 
system.  The  phases  are  modeled  according  to  the  phase  stabilities,  measured  thermodynamic 
properties,  and  the  measured  transus  points  (including  liquidus,  solidus,  eutectic,  etc).  In  our 
research,  differential  scanning  calorimetry  (DSC)  is  used  to  establish  thermal  events  associated 
with  phase  transformations.  X-ray  and  electron  diffraction  determine  structure  and  analytical 
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electron  microscopy  may  all  be  used  to  validate  equilibrium  phase  boundaries  calculated  with  the 
CALPHAD  techniques. 

For  alloy  systems,  phases  may  be  formally  treated  using  a  sub-lattice  concept.  For 
ternary  non-ideal  solutions,  their  Gibbs  energy,  using  a  sub-regular  model,  may  be  expressed  in  a 
Redlich-Kirst-Muggianu’s  formula. 


RTp,  tox,  +  ix,x,  I  ‘i,./‘’(x,  -  x;‘>)  ^  |x,x,x,  ^  ♦  £.n  (2) 


The  specific  models  employed  and  the  interaction  parameters  L  between  elements,  e.g.,  i  and  j 
for  each  phase,  are  compiled  in  a  database.  Retrieving  the  data  from  the  database,  a  software- 
package  Thermo-Calc  can  then  be  used  to  determine  the  stable  phases  present  in  a  system  as  a 
function  of  composition,  temperature  and  pressure.  The  stable  phases  are  determined  from  the 
relative  minimum  in  the  Gibbs  free  energy  term  for  each  phase  possible  in  the  system. 

Regarding  the  Al-Cu-Mg-(Ag)-X  (X=Si,  Zn,  Ce...)  system,  well-developed  databases  for 
Al-Cu-Mg,  Al-Cu-Si,  Al-Mg-Si,  Al-Cu-Mg-Zn  etc.  are  available.  The  elements  including  Al,  Si, 
Mg,  Zn,  Ce  etc.  are  considered  in  a  light-alloy  database  called  “COST507”  [13]  developed  by  a 
European  group  with  participants  of  14  EU  countries.  We  also  have  been  using  the  database  TT- 
A1  [14]  that  is  useful  for  solidification  studies  near  the  Al-rich  comer  of  the  system.  Effort  was 
also  made  for  the  construction  of  an  Al-Cu-Mg-Ag  database  [15].  However,  the  most  significant 
information  for  our  alloy  design  purpose,  particularly  in  the  Al-rich  comer,  were  either  missing 
or  incomplete; 

1)  Both  Q  and  6',  the  selected  strengthening  phases,  have  not  been  modeled; 

2)  The  S  phase  was  modeled  but  Ag-partitioning,  the  controlling  factor  for  phase 

competition  with  Q  and  0',  was  not  incorporated. 

3)  The  effects  of  the  trace  element,  e.g..  Si,  Zn,  Ce  etc.,  have  not  been  characterized. 

The  first  critical  component  needed  to  optimize  this  class  of  alloys  is  to  control  the  bulk 
composition  so  that  it  lies  within  the  (a-FCC  +  0'  +  Q)  phase  field  in  order  to  avoid  S-phase 
precipitation  during  aging  treatment  (since  S-particles  are  less  effective  in  strengthening  and  act 
to  remove  Mg  and  Ag  from  the  Q  interface  which  leads  to  the  dissolution  of  the  Q  plates).  To 
optimize  the  Al-Cu-Mg-Ag  quaternary  system  the  (a+0'+Q)/(a+0'+Q+S)  phase  boundary  must 
be  accurately  located  and  determined  as  a  function  of  bulk  composition,  artificial  aging  time  and 
temperature.  The  databases  can  also  be  used,  as  shown  below,  for  evaluation  of  some  properties 
that  are  reflected  in  the  thermodynamics  of  the  system. 

Formation  of  the  Desired  Precipitate  Structure;  In  the  Al-Cu-Mg-Ag-X  system,  0'- 
precipitates  evolve  from  clustering  of  Cu  on  {100}„-planes  in  Al-FCC.  The  nucleation  of  Si- 
plates  is  assisted  by  clustering  or  short-range  ordering  of  Cu-Mg  and/or  Ag-Mg,  that  act  as 
precursors,  as  observed  in  [9,10].  Since  the  formation  of  the  undesired  S-phase  will  consume 
Cu,  Mg  and  Ag  from  the  matrix,  this  phase  is  deleterious  for  the  nucleation  of  the  0'-  and 
particularly  Q-phase  and,  hence,  for  the  formation  of  the  desired  (0^+Si)  structure.  Dislocations 
introduced  by  pre-age  plastic  transformation,  the  T8-treatment,  were  also  found  to  favor  the 
nucleation  of  0'  [16].  The  strategy  for  our  plan  for  increasing  the  nucleation  rate  of  the  Q  phase 
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focused  on  the  enhancement  of  the  formation  of  precursors  (clusters  or  SROs)  by  trace  additions. 
The  desired  (O’+Q)  structure  with  a  fine  and  uniform  distribution  will  eventually  be  produced 
through  rigid  control  of  Cu  and  Mg  contents,  trace  element  additions  and  using  refined  thermo¬ 
mechanical  processing. 

Thermal  Stability  of  the  Desired  Precipitate  Structures:  The  coarsening  behavior  of  the  Q 
plates  in  Al-Cu-Mg-Ag  has  been  extensively  studied  in  our  lab  [11].  The  coarsening  or 
thickening  rate  is  slower  when  the  thickness  of  Q-plates  is  less  than  ~  ISnm.  This  may  be 
attributed  to  the  inability  to  form  the  misfit  compensating  in-situ  dislocations  at  the  edge 
interface  between  the  plates  and  the  matrix.  The  vacancy  (dilation)  strain  associated  with  the 
(-9%)  misfit  between  the  single  _  unit  cell  high  £2-plates  and  2  {1 1 1}„  plane  spacing  makes  the 
nucleation  of  thickening  ledges  difficult,  particularly  during  the  coarsening  stage  where  the 
driving  forces  will  be  much  lower  than  during  growth.  When  the  plates  are  thicker  than  the 
critical  value,  acceleration  in  thickening  occurs  due  to  the  relaxation  of  the  strain  field  with  the 
formation  of  the  misfit  dislocations.  Another  important  finding  is  the  trace  element  segregation 
in  the  broad  interfaces  of  Q  plates. 


More  generally,  there  are  two  issues  for  the  thermal  stability  of  a  precipitate  structure 
containing  more  than  one  2"**  phase  -  phase  coarsening  and  phase  competition  from  other  phases. 
Coarsening  of  a  single  2"'*  phase  is  thermodynamically  attributed  to  the  tendency  for  decreasing 
the  total  interfacial  energy  between  the  2"^  phase  and  the  Al-matrix.  The  interfacial  energy  may 
consist  of  the  chemical  energy,  as  described  by  the  bond-breaking  model,  and  the  associated 
strain  energy.  Trace  element  additions  may  have  two  possibly  “useful”  ways  of  partitioning: 
either  segregating  at  the  interfaces  or  being  soluble  in  the  2"**  phase.  The  former  alters  the 
chemistry  of  the  interfaces,  e.g.  Ag  layers  observed  by  Z-contrast  HRTEM  [1 1],  possibly  leading 
to  a  lower  interfacial  energy.  The  partitioning  into  tihe  2”‘'  phase,  like  the  Zn-partitioning  in  Ti 
phase  [17],  changes  the  lattice  constants  and  may  hence  modify  the  misfit  strain.  It  was  also 
shown  that  Ag  segregates  to  S  phase  in  an  Al-Cu-Mg-Ag  alloy  [10].  On  the  other  hand,  different 
from  the  simple  alloys  in  the  two-phase  field,  the  coarsening  behavior  of  higher-order  alloys  like 
the  Al-Mg-Cu-(Ag)-(Si)-X  system  containing  more  than  one  2"**  phase  (like  Q  ,  0'  and  S)  is  also 
affected  by  the  competition  from  other  2"**  phases  for  the  solute  Cu,  Mg  and/or  Ag  that  they 
share.  For  example  at  service  temperatures  (<150°C)  usually  lower  than  aging-treatment 
temperatures  >200®C,  the  driving  force  for  formation  of  S-phase  will  inevitably  become 
significant.  The  formation  of  S-phase  consuming  the  solute  Mg  and  Ag,  that  are  stabilizers  for  Q, 
and  Cu  in  the  A1  matrix  and,  hence,  make  both  0'  and  Q  particles  less  stable  leading  to  potential 
dissolution.  Therefore,  our  objective  here  is  not  only  to  lower  the  coarsening  rate  of  the  Q  and  0'- 
plates  but  also  to  prevent  the  S-phase  from  jeopardizing  the  stability  of  the  Q  plates. 

First  Principle  Calculations:  Crystal  Structure,  Elastic  Constants  and  Interfacial  Enei^es: 
Although  the  precipitate  phases  (Q,  0',  S,  etc.)  in  this  system  have  been  known  for  years,  some 
fimdamental  properties  such  as  crystal  structure,  enthalpy  of  formation,  deformation  behavior 
and  interfacial  energy  between  these  phases  and  the  A1  matrix  are  not  available  because  of 
difficulties  in  experimental  measurement.  These  data  are  necessary  for  any  reliable  quantitative 
thermodynamic  and  kinetic  analysis  of  this  system. 
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First  principle  total  energy  calculations  are  being  employed  for  this  purpose.  We  have 
been  using  the  Vienna  ab  initio  simulation  program  (VASP)  [18]  where  the  density  functional 
theory  plane-wave  pseudo-potential  method  is  applied.  Mainly  LDA  (Local-Density 
Approximation)  is  used  for  electronic  exchange  correlation  potentials.  GGA  (general  gradient 
approximation),  such  as  PWE  and  PW91  functionals,  are  also  being  tested.  The  ultrasoft  pseudo¬ 
potentials,  supplied  with  the  VASP-package,  are  utilized.  The  GGA  did  not  produce  any 
significant  changes  in  niost  of  our  results.  The  basis-set  cut-offs  were  set  to  “high  accuracy”  and 
k-point  meshes  are  set  by  “Monkhorst-Pack”  (or  T-centered)  setting  from  8x8x8  to  32x32x32 
(or  otherwise  indicated)  depending  on  the  size  of  the  unit  cell  concerned. 

Crystal  Structures  and  Formation  Enthalpy:  The  first  principle  energy  calculations  are  used 
to  validate  the  proposed  crystal  structures  and  to  calculate  the  formation  energies  of  various 
precipitates.  It  offers  some  important  data  that  are  unavailable  from  an  empirical  approach.  The 
formation  enthalpies  obtained  may  be  used  for  input  in  the  CALPHAD  database  development. 
For  a  binary  phase  ApB^,  the  formation  enthalpy  is  defined  as  the  energy  gain  or  loss  with  respect 
to  the  bulk  constituents  A  and  B  at  their  equilibrium  lattice  constants 


where  and  E“KA)and  E“KB)  are  the  energies  of  the  second  phase  ApB,  and  the 

constituents  A  and  B.  x^  =  p/(p+q)  and  Xg  =  q/(p+q)  are  the  concentrations  of  A  and  B.  The 
formation  energy  corresponds  to  geometrically  fully  relaxed  compounds,  i.e.,  the  structures  are 
optimized  (consistent  with  the  symmetry  of  the  structure)  with  respect  to  unit  cell  vectors,  atomic 
displacements  and  volume  of  the  unit  cell. 

0'  and  0-phases:  The  tetragonal  crystal  structures,  including  space  group,  internal  atomic 
coordinates  and  lattice  constants  of  0'  and  0-phases  have  been  known  for  some  time.  Our 
calculations  gave  the  formation  enthalpies  and  validated  the  lattice  constants  as  shown  in  Table  I. 

Table  I.  Experimental  and  Calculated  Lattice  Constants  and  Calculated  Enthalpy  of  0'  and  0- 
phases 


Phase 

Lattice  Constants- 
Experimental,  A 

Lattice  Constants- 
Calculated,  A 

Enthalpy  of  Formation- 
Calculated,  eV/atom 

0' 

A=  4.04,  c=  5.80 

4.0756,4.0756,5.8294 

-0.213 

0 

A=6.063,  c=4.87 

6.074,6.074,4.850 

-0.168 

Q  Phase:  Several  crystal  structures  for  the  Q  phase  have  been  proposed.  Auld  suggested 
a  monoclinic  structure  [4].  A  widely  accepted  structure  [7]  is  orthorhombic  with  space  group 
(Fnrnim).  Based  on  analysis  of  the  convergent-beam  electron  diffraction  (CBED)  patterns,  Garg 
and  Howe  [8]  advocated  a  tetragonal  structure  (I4/mcm)  that  is  only  a  metastable  {lll}a  variant 
of  the  equilibrium  0  phase.  Table  II  lists  our  calculated  results  using  VASP.  It  shows  that  the 
proposed  monoclinic  structure  has  a  positive  formation  enthalpy  and  hence  is  unreasonable.  Both 
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the  tetragonal  and  the  orthorhombic  structures  have  formation  enthalpies  that  are  close  to  that  of 
the  equilibrium  6  phase.  According  to  the  COST  507  database  that  was  compiled  via  the 
CALPHAD  method,  the  formation  enthalpy  of  the  equilibrium  0  phase  is  -161  meV/atom. 
Experiments  have  yielded  values  from  -245  to  -135  meV/atom  [19].  Due  to  the  similarity  of  the 
calculated  formation  enthalpy  of  Q  to  that  of  equilibrium  0,  Q  may  be  regarded  as  a  metastable 
{11 1}„  variant  of  0. 


Table  II.  VASP  Calculated  Formation  Enthalpies  and  the  Vector  Lengths  (unit  cell)  of  Proposed 
Crystal  Structures  for  the  Q  phase. 


Proposed 

H,  eV/atom 

Vector  Length,  A 

Equilibrium  0 

Tetragonal  (I4/mcm): 
6.063,4.87 

-0.168 

6.074,6.074,4.850 

Auld’s 

Monoclinic: 
a  =  b:4.%,c:8.48,  120 

0.195 

4.930,4.930,9.122 

Knowles’ 

Orthorhombic 

(Fmmm);4.96,8.59,8.46 

-0.167 

4.8946,8.7120,8.4151 

Garg  &  Howe’s 

Tetragonal  (I4/mcm): 
6.066,4.% 

-0.169 

6.073, 6.073,4.855 

$-phase  (AljCuMg):  There  are  also  several  models  for  the  S  phase  based  on 
experimental  observations.  The  widely  accepted  structure  is  base-centered  orthorhombic  with 
space  group  Ctncm  and  lattice  constants  a=4.03,  b=9.30  and  c=7.08.  However,  there  is  some 
dispute  as  to  internal  atomic  coordinates.  Our  calculation,  verifying  a  similar  study[191,  showed 
that  the  original  model  by  Perlitz  and  Westgrin  leads  to  a  formation  enthalpy  of  -0.210  eV/atom 
while  exchange  of  the  Cu  and  Mg  in  internal  atomic  coordinates,  as  suggested  by  Radmilovic  et 
al  [20],  yielded  a  positive  formation  enthalpy  of  +0.289  eV/atom.  By  neglecting  Ae  entropy 
contribution  (mostly  from  lattice  vibration),  this  indicates  that  the  Perlitz  &  Westgrin  model  is 
reasonable. 

Interfacial  Energy  Between  the  2nd  Phases  with  AI  Matrix:  Interfacial  energy  of  a  second 
phase  in  the  matrix  is  the  key  parameter  for  stability  of  precipitates  and  hence  aid  in  determining 
the  creep  resistance  of  the  alloys.  Using  the  First  Principle  approach,  the  coherent  interfacial 
energy  between  plate-like  precipitates  and  the  matrix  can  be  calculated  by  setting  up  an 
appropriate  super-cell  containing  the  interfaces.  The  interfacial  energy  Yi  may  be  defined  as 

V  B  where  E„  and  denote  the  formation  energy  of  the  precipitates  and 

>>  A  a  p 

^(uv) 

the  matrix  per  atom,  the  n„and  n^  the  number  of  the  atoms  contained  in  the  supercell,  and  A(„v)  is 
the  area  of  the  interface  of  the  supercell.  Figure  2  shows  the  supercell  construction  for 
calculation  of  the  S2(0m) — Al  (fee)  interfacial  energy  parallel  to  the  Q-habit  plane  (lll)a  where 
the  coinciding-site  lattice  in  the  interface  was  realized. 
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[121]„  //[OOlJej, 


[111]„//[110U 
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Figure  2.  Supercell  construction  for  calculation  of  the  interfacial  energy  between  Q  (0m)  and  A1  (111).  The 
top  is  the  projection  on  YZ-plane  and  the  middle  is  the  projection  on  XZ-plane.  The  bottom  shows  the 
orientation  match  between  the  two  phases.  The  pink  balls  represent  Al-atoms  and  the  brown  Cu-atoms  in  the 
Q-phase,  and  the  blue  are  atoms  in  the  Al-(fcc). 


The  ideal  calculation  is  supposed  to  converge  with  the  k-points  meshing,  energy  cut-off 
and  size  of  the  supercell.  A  satisfactory  result  of  170  mJ/m^,  has  been  obtained  for  the 
0' — Al(fcc)  broad  interface  parallel  with  the  habit  plane  {100}.  The  result  of  another 
independent  calculation  by  V.  Vaithyananthan  et  al.  was  235  mJ/m^  [21].  The  calculations  for 
other  second  phases  are  currently  running  on  our  computer.  It  is  noted  that  since  the  internal 
atomic  coordinates  close  to  the  interfaces  are  relaxed,  the  interfacial  energy  obtained  consists  of 
contributions  from  both  chemical  “bond-  breaking”  effect  and  the  misfit  strain  energy.  Again  the 
results  exclude  the  change  in  the  entropy  due  to  the  interfaces.  Interfacial  energy  calculations 
will  also  be  used  to  select  trace  alloy  additions  that  may  reduce  the  interfacial  energy  and 
increase  the  thermal  stability  of  the  strengthening  precipitates. 

Deformation  Behavior  of  Clusters  and  Precipitates:  Quantitative  data  of  the  elastic  constants 
for  the  2nd  phases  are  significant  and  fundamental  for  determining  the  strengthening 
mechanisms.  These  parameters  also  play  a  role  in  determining  the  habit-plane  and  the  evolution 
of  equilibrium  shapes  of  the  clusters  and  the  2nd  phase  particles  since  they  determine  the  elastic 
strain  energy  associated  with  lattice  misfit. 

Direct  measurement  of  the  2“*  phase  elastic  constants  of  second  phases  is  usually 
impractical.  Quantitative  analysis  based  on  the  elastic  constants  typically  employ  very  rough 
estimates.  First  principle  calculations  can  be  employed  to  calculate  the  total  energy  of  a  crystal 
structure  at  0  K.  The  2“*  order  elastic  constants  of  the  crystalline  solid  at  0  K  may  be  defined  and 
calculated  according  to  the  variation  in  the  energy  from  its  optimized  lattice  structure  as 
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C  «  ( — Zr — )  with  i,i.k,l  =  1,2,  or  3  where  E  denotes  the  total  energy,  Ejj  the  strain 
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tensor.  Using  “engineering”  notations,  Hooke’s  law  is  applied;  a,  *  ^  where  Oy  is  the 

stress  tensor  and  Sy  the  elastic  compliance  tensor.  By  applying  strains  in  certain  patterns,  we  can 
then  retrieve  the  variation  of  the  energy  or  the  resulting  stress  tensor  and  calculate  the  elastic 
constants  Gyu  or  Syu.  Table  III  lists  the  calculated  results  for  Al,  Cu,  Ag,  and  Mg  (hep).  We  note 
that  the  results,  although  strictly  for  0  K,  are  quite  consistent  with  experimental  measurements. 
The  calculated  results  for  0’,  0,  Q  and  S-phase  are  listed  in  Table  IV . 


Table  III.  Elastic  constants  calculated  using  VASP  and  obtained  from  experiments  at  room  temperature 


for  some  element  crystals 

1  Elastic  Constants,  10*0  N/m2 

Calculated 

Experiments 

Al  fee  C„=9.30,  C, 2=6.38,  C44=3.1 

J,  B=7.35  C„=10.8  Cif=6.13  C44=2.85,B=7.2 

Cu  fee  C„=17.9,C,2=12.4,C44=8.03 

B=14.27  C„=16.84,C, 2=12.14, C44=7.54,B=13.7 

Ag  fee  C„=10.93,C,2=8.39,C44=3.7 

:,  B=9.24  C„=12.4,C,2=9.34,C44=4.61,B=10.36 

Mg  hep  C„=5.39,C33=4.78,C,2=2.2 

C^=1.75 

9,  C, 3=3.26,  C„=5.64,  C33=5.81,  Ci2=2.32,  C, 3=1.81 
C44=1.68 

Table  IV.  Elastic  constants  calculated  using  VASP  for  some  second  phases 


Phase 

0’ 

Q/0/0„ 

S 


Structure  Elastic  Constants,  10^^  N/m^ 

tetr.  Cii=  19.24  ,  C33  =  15.62,  €,2=2.91,  €,3  =  6.54,  €44=8.35,  Cgs  =  4,36 

tetr.  Cii=18.72,  5.59,  =  5.36  €33=18.30,  C44  =  3.68,  6.27 

orth.  C„=14.5,C22=15.3,C33=14.3,Ci2=  2.95,Ci3=6.60,C44=5.31,C5;p7.29,Cs,p3.29 


Alloy  Partitionii^  and  the  Al-Cu-Mg-(Ag)  Phase  Di^am:  A  series  of  isothermal  diagrams 
for  the  Al-Cu-Mg  ternary  system  are  shown  in  Figure  3.  It  is  important  to  note  that  the  (a  + 
0)/(a  +  0  +  S)  phase  boundary  shifts  to  higher  Mg  contents  with  increasing  temperature,  which  is 
found  at  approximately  0.15  "Vo  and  0.50  "*/o  for  200°C  and  300‘’C,  respectively.  This  indicates  a 
higher  solid  solubility  for  Mg,  which  limits  the  driving  force  for  the  formation  of  the  S-phase,  as 
temperature  increases.  A  compromise  must  be  reached  concerning  the  Mg  levels  in  the  proposed 
alloys.  In  order  to  avoid  S-phase  precipitation  the  Mg  contents  must  be  kept  to  very  low 
amounts  (~0.10'^/o)  and  precisely  controlled  in  these  alloys.  However,  the  allowable  Mg 
content  should  be  maximized  in  order  to  optimize  the  precipitation  of  £2,  recalling  that  it  is 
necessary  for  nucleation  at  moderate  temperatures.  As  a  point  of  interest  current  experimental 
Al-Cu-Mg-Ag-X  alloys  lie  within  the  (a  +  0  +  S)  phase  field  for  artificial  aging  at  moderate 
temperatures  (150-250“C)  according  to  the  calculated  ternary  isothermal  sections.  Also,  it  is 
important  to  note  that  this  diagram  does  not  account  for  the  Q  phase.  Rather,  Q  precipitation  is 
treated  as  a  {1 1 1}„  variant  of  equilibrium  0,  as  previously  mentioned. 
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Figure  3.  Calculated  montage  of  isothermal  sections  of  the  Al-Cu-Mg  ternary  system. 

We  have  verified  Ringer  et  al’  discovery  [10]  of  Ag  partitioning  to  the  S-phase 
precipitates,  as  shown  in  Fig.  4.  Systematic  quantitative  analysis  of  the  Ag  partitioning  to  S- 
phase  precipitates  is  in  progress.  This  analysis  will  provide  the  chemical  composition  of  S-phase 
precipitates  as  a  function  of  artificial  aging  temperature  so  that  the  energy  descriptions  for  this 
phase  in  the  Thermo-Calc  database  will  account  for  the  trace  additions  of  Ag,  rather  than  the 
stoichiometric  ternary  compound  as  modeled  in  the  current  database  [15].  The  quantitative 
information  could  also  aid  in  the  calculation  of  the  interfacial  energy  between  S-phase 
precipitates  and  the  Al-rich  matrix.  Due  to  the  complexity  of  the  interface  structure,  first 
principle  calculation  of  the  interfacial  energy  is  impractical. 
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Figure  4.  a)  Representative  EDS  spectra  collected  from  within  a  S-phase  precipitate  aged  at  250“C  for 
lOOOh.  b)  High-energy  regime  of  EDS  spectra  illustrating  the  partitioning  of  Ag  to  the  S-phase 
precipitate  relative  to  the  surrounding  matrix. 

Effect  of  Silicon  Content  on  Precipitation  of  the  Q  Phase:  Despite  the  rigorous  studies 
performed  on  Al-Cu-Mg  alloys  from  1929-1954,  the  Q  phase,  in  the  absence  of  Ag,  was  not 
identified  until  1990.  Revisiting  a  summary  of  the  studies  performed  on  Al-Cu-Mg  alloys  from 
1929-1954  confirms  that  the  Si  content  of  the  majority  of  the  alloys  studied  exceeded  O.l’^/o*  For 
many  of  these  alloys.  Si  was  intentionally  added  to  increase  strength  (arising  from  enhanced 
matrix  precipitation  of  the  strengthening  phases).  These  Si  levels  may  have  led  to  the 
suppression  of  Q  precipitation  in  these  early  Al-Cu-Mg  alloys,  and  Q  was  not  observed  until 
microstructural  characterization  was  performed  on  high  purity  alloys. 

Our  experimental  data  [22]  has  demonstrated  that  Q  precipitation  is  suppressed  in  alloys 
with  Mg/Si  atomic  fraction  ratios  of  s  2.0,  regardless  of  Si  content  and  indicates  that  Q 
precipitation  is  not  dictated  by  a  critical  Si  content,  but  rather  by  the  Mg/Si  ratio  of  the  alloy.  It  is 
important  to  note  that  the  critical  Mg/Si  ratio  of  2.0  is  consistent  with  the  expected  bonding 
between  these  two  elements  based  upon  atomic  valency.  Similarly,  the  MgjSi  (P')  precipitate 
formed  in  Al-Mg-Si  alloys  has  an  atomic  Mg/Si  ratio  of  2.0. 

Our  results  demonstrate  that  trace  amounts  of  Si  were  responsible  for  limiting 
precipitation  in  low  Mg-containing  alloys.  They  indicate  that  an  atomic  Mg/Si  ratio  of 
approximately  2.0  must  be  overcome  for  Q  nucleation.  This  sensitivity  of  Q  precipitation  to  Si 
content  may  explain  the  absence  of  Q  in  early  work  on  Al-Cu-Mg  alloys  and  must  be  taken  into 
consideration  in  our  design  of  Al-Cu-Mg-(Ag)  alloys  for  high  strength  and  creep  resistance  at 
moderate  temperatures. 

Precursors  and  Non-Precursors  for  the  2"^  Phases  —  SRO  of  Solutes:  To  achieve  our 
objective  to  produce  ‘Tine,  uniform  precipitate-structures”,  the  strategy  is  to  focus  on  boosting 
the  nucleation  rates  of  the  favored  strengthening  precipitates  through  enhancement  of  the 
formation  of  their  precursors  by  trace  addition  and/or  elimination  of  harmful  impurities. 


Figure  5.  The  short-range  ordering  parameters  (a)  of  Mg-Si,  and  (b)  a^’®'  of  Al-Si  (Warren- 
Cowley  parameter  /  X/^b  where  ^^’®  denotes  the  pairing  probability  between  A  and  fi¬ 
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atoms)  in  Al-Mg-Si  fcc-Al  phase  calculated  using  quasi-chemical  approximation  and  CALPHAD 
databases  -  COST2. 

The  negative  effect  of  Si  in  Al-Cu-Mg  alloys  mentioned  above  can  be  ascribed  to  the 
formation  of  short-range-ordering  (SRO)  like  Mg-Si,  that  we  call  “non-precursors”,  in  the  Al- 
FCC  solid  solution.  It  consumes  the  solutes  needed  for  the  precursors  for  the  Q-phase  and  hence 
jeopardizes  its  nucleation.  We  have  set-up  a  method  to  evaluate  SRO  in  order  to  rule  out 
potential  deleterious  elements  like  Si  in  Al-Cu-Mg.  The  method  is  based  on  the  quasi-chemical 
approximation  with  consideration  of  the  1**  neighbor  interaction  [23].  The  interaction  parameters 
between  each  pair  of  elements  are  calculated  using  the  heat  of  mixing  of  the  alloy  system  under 
CALPHAD  spirit.  Using  a  sub-regular  solution  model  as  shown  in  Eq.  (1)  for  a  ternary  FCC 
solution  Al-Mg-Si,  the  heat  of  mixing  may  be  written  as 


= 


=2  J 

M  j  y~..i 


0''\j -x^r  ^x,x,x,'^0\2.3 


d  1*2^ 
T — 

dr 


(4) 


where  the  interaction  L-parameters  can  be  retrieved  from  available  Thermo-Calc  databases.  Fig. 
5  shows  the  iso-values  of  the  SRO  parameters  of  Mg-Si  and  Al-Si  in  ternary  system  Al-Mg-Si  - 
FCC  solutions  at  200"C.  As  indicated  by  very  negative  Warren-Cowley  parameters  when 
compared  to  the  mostly  positive  values  for  Al-Si  in  the  composition  range  interested.  Si  has  a 
strong  tendency  to  form  Mg-Si  SRO  in  the  FCC  (a)  phase  and  hence  to  jeopardize  the  formation 
of  £2. 


Conversely  some  trace  element  additions  may  have  a  very  positive  effect  on  nucleation  of 
the  desired  {Q  and  0'}-phases  as  shown  above.  In  the  Al-Cu  binary  alloy,  Cu  tends  to  cluster  on 
{100}-planes  in  Al-Cu  solutions.  These  clusters  evolve  into  0"  and  0'  with  {100}a  habit  planes. 
However,  the  addition  of  Mg  with  Cu:Mg=10:l  leads  to  the  tendency  to  form  Q  plates  on  {111}. 
As  indicated  in  the  3D-atomic  probe  (3D-AP)  and  HRTEM  observation  [9,10],  this  can  be 
attributed  to  the  clustering  or  short-range  ordering  (atom  aggregates)  as  the  precursors  formed 
prior  to  the  precipitation  of  the  Q  phase.  Analysis  of  elastic  strain  energy  caused  by  these  plate¬ 
like  clusters  may  elucidate  what  Mg  and/or  Ag  clusters  may  function  as  the  precursors  for  £2 
{lll}-phases.  Assuming  dilute  solute  approximation,  we  have  set-up  a  method  to  calculate  the 


elastic  strain  energy  dependence  on  the  orientation  <hkl>a  of  the  plate-like 
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clusters  embedded  in  the  Al-FCC  matrix.  The  elastic  constants  Cy^  were  obtained  using  first- 
principle  calculations.  Fig.  6  shows  the  parametric  plots  of  the  orientation  dependence  of  elastic 
strain  energy  for  some  clusters.  It  shows  that  the  Cu-FCC  (Cii=17.90,  Ci2=12.46,  €44=8.03) 
clusters  has  the  highest  strain  energy  as  they  are  parallel  with  {lll}a,  and  the  lowest  energy 
strain  energy  when  parallel  with  {100}„.  This  explains  the  formation  of  Cu  plate-like  clusters  on 
the  planes  {100}„,  as  experimentally  observed  Cu  GP-zones  in  Al.  The  hypothesized  Mg-FCC 
structure  (Cij=4.76,  Ci2=  3.09,  €44=2.27)  is  similar  to  Cu.  Although  the  possible  effect  of 
interfacial  energy  on  the  shape  of  Mg  clusters  may  explain  why  Mg-GP  zones  are  spherical,  the 
calculation  indicates  that  pure  Mg  atoms  are  difficult  to  form  plate  clusters  on  {1 1 1}„  where  they 
have  the  largest  strain  energy.  Fig.  6  also  shows  a  preliminary  result  for  a  “quasi-disordered” 
FCC-mixture  of  Cu  and  Mg  (3:1).  It  shows  that  this  mixture  has  the  lowest  strain  energy  when 
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parallel  with  .  This  indicates  that  it  is  the  mixture  clusters  of  Cu-Mg,  that  most  likely 

form  on  {lll}o„  that  functions  as  precursors  for  £2  precipitates.  The  strong  calculated  pairing 
probability  or  negative  SRO-parameter  for  Mg-Cu  supports  this  point 


Figure  6.  parametric  plot  of  the  orientation  {M/}-dependence  of  elastic  strain  energy  (under 
dilute  approximation)  caused  by  embedding  a  plate-like  FCC  cluster  in  Al-matrix  oriented 
with  respect  to  the  Al-cubic  coordinates,  a)  Cu-FCC;  b)  Mg-FCC  (postulated);  and  c)  a 
Cu3Mgl-FCC  mixture.  Insert  d)  indicates  the  cubic  coordinates  of  the  FCC-matrix. 

We  have  established  two  factors  to  identify  potentially  “useful”  trace  elements  for 
forming  precursors  and/or  “harmful”  impurities  that  form  non-precursors:  the  SRO  parameter 
that  gauges  probability  for  potential  elements  to  form  either  SROs  or  clusters  (atom  aggregates), 
and  the  orientation-dependence  of  the  strain  energy  that  controls  the  habit-plane  selection  of 
these  plate-like  aggregates.  Further  SRO  analysis  depends  on  the  more  complete  CALPHAD 
databases  we  are  currently  developing  for  Al-Cu-Mg-Ag-(X).  The  calculation  of  the  elastic  strain 
energy  for  the  mixture  or  alloy  clusters  has  employed  the  elastic  constants  obtained  by  First 
Principle  calculations.  Size  limitation  of  the  (“quasi-disordered”)  supercell  constructed  due  to 
computer  capability  does  not  allow  for  more  realistic  modeling  of  the  disordered  mixtures.  Very 
recently,  a  robust  technique  was  developed  for  calculation  of  elastic  constants  of  disordered 
mixtures  or  alloys  [24,25].  It  is  based  on  quasi-chemical  approximation  with  consideration  of 
atomic  pair-interaction  with  the  1**  and  2“*  neighbors.  The  elastic  constants  are  calculated  from 
the  axis-synunetric  force  constants  of  a  virtual  pair  potential  constructed  for  atoms  interaction  in 
the  disordered  alloys.  The  virtual  potentials  can  be  evaluated  using  curve-fitted  data  for  pair 
bond-lengths  from  measured  mean  values  or  lattice  parameters  of  the  alloy  crystals  and/or  from 
the  known  elastic  constants  of  crystals  (with  the  same  structure)  of  the  terminal  pure  elements. 
The  technique  can  be  applied  to  binary  and  ternary  cubic  alloys.  It  is  noted  that  very  often  when 
terminal  single  crystals  with  the  same  structure  are  unavailable,  the  first  principle  calculation  can 
be  used  for  that  input.  Therefore,  our  future  work  will  focus  on  combining  the  first  principle 
calculations  with  this  “virtual  potential”  technique  for  the  orientation  selection  analysis  of  the 
embedded  mixture-clusters  e.g.,  Ag-Cu,  Ag-Mg,  Mg-Ag-Cu,  Zn-Cu  etc,  that  are  critical  as 
potential  precursors  for  the  desired  £2  and/or  0'-phases. 


The  Effect  of  Cold  work  on  the  Precipitation  and  Coarsening  Behavior  of  Q  and  0':  The 
effect  of  plastic  deformation  on  the  microstructural  evolution  of  an  Al-5.0Cu-0.5Mg  alloy  was 
investigated  in  order  to  study  how  plastic  deformation  affects  the  competitive  nucleation,  growth 
and  coarsening  behavior  of  Q  and  0'  in  the  absence  of  Ag  [26].  Hardness  measurements  and 
quantitative  precipitate  analysis  were  performed  on  specimens  that  were  water  quenched  from  a 
solution  heat  treatment,  stretched  either  0%  or  6%  and  immediately  aged  at  ambient  temperature 
or  artificially  aged  at  200  °C  or  250  °C  for  times  up  to  3000h.  Quantitative  transmission  electron 
microscopy  (TEM)  was  used  to  characterize  Q  and  6'  number  density,  diameter  and  thickness  as 
a  function  of  pre-age  mechanical  stretch  and  artificial  aging  condition.  Age  hardening  curves  for 
naturally  and  artificially  aged  specimens  revealed  an  increase  in  hardness  corresponding  with  a 
pre-age  stretch.  Quantitative  TEM  verified  an  increase  in  the  number  density  and  a  refinement 
of  precipitates  for  both  Q  and  0'  between  the  0%  and  6%  stretch  condition  for  those  samples 
artificially  aged.  When  aged  at  200  "C,  0'  exhibited  superior  coarsening  resistance  relative  to  the 
Q  phase.  The  quantified  Q  coarsening  kinetics  was  greater  than  similar  Ag-containing  alloys 
indicating  that  Ag  may  change  the  interfacial  energy  or  the  strain  energy  of  the  Q  precipitates. 
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